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While the Sauter-Schwinger effect describes nonperturbative electron-positron pair creation from 
vacuum by a strong and slowly varying electric field -Batrong via tunneling, the dynamically assisted 
Sauter-Schwinger effect corresponds to a strong (exponential) enhancement of the pair-creation 
probability by an additional weak and fast electric or electromagnetic pulse Bweak- Using the WKB 
and worldline instanton method, we find that this enhancement mechanism strongly depends on 
the shape of the fast pulse. For the Sauter profile 1 / cosh^ (cut) considered previously, the threshold 
frequency cUcrit (where the enhancement mechanism sets in) is basically independent of the magnitude 
Bweak of the Weak pulse—whereas for a Gaussian pulse exp(—cu^t^), an oscillating profile cos(cut) or 
a standing wave cos(cut) cos(fca:), the value of cucrit does depend (logarithmically) on Bweak/Batrong- 

PACS numbers: 12.20.-m, ll.15.Tk 


I. INTRODUCTION 

In contrast to the perturbative realm of quantum 
field theory, where we possess well-established methods 
for calculating observables such as scattering cross sec¬ 
tions, our understanding of nonperturbative effects is still 
rather incomplete. In quantum electrodynamics (QED), 
a prominent example for such a nonperturbative phe¬ 
nomenon is the Sauter-Schwinger effect m correspond¬ 
ing to the creation of electron-positron pairs out of the 
quantum vacuum by tunneling induced by a strong elec¬ 
tric field E. For constant electric fields E, the leading- 
order pair-creation probability scales as 



where ig is the charge and m the mass of the elec¬ 
trons and positrons. Since the critical field strength 
Es ~ 10^®V/m is extremely large, this prediction has 
not been conclusively experimentally verified yet |5] . 

Several years ago, it was found [7] that the above pair- 
creation probability P^+e- can be drastically enhanced 
by superimposing the constant (or slowly varying) strong 
field E with a weaker time-dependent pulse—the dynam¬ 
ically assisted Sauter-Schwinger effect; see also Refs. [8l- 
Ml- Interestingly, this enhancement mechanism is al¬ 
ready operative for frequency scales co far below the mass 
gap 2mc^ separating the Dirac sea from the positive con¬ 
tinuum. In the following, we study the dependence of 
the dynamically assisted Sauter-Schwinger effect on the 
shape of the additional time-dependent pulse. To this 
end, we employ the WKB and the worldline instanton 
methods in order to compare a Gaussian pulse and an 
oscillating or standing-wave profile with the Sauter pulse 
considered previously. 
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II. RICCATI EQUATION 

Let us briefly review the main steps of the derivation 
(see e.g. Ref. [15] for a more detailed introduction). For 
simplicity, we start with the Dirac equation in 1+1 di¬ 
mensions {h = c = 1) 

(iY + iqA^] + to ) • -0 = 0 , ( 2 ) 

with the two-component spinor i/). The purely time- 
dependent electric field E{t) is treated as an external 
background field and can be described by the vector po¬ 
tential = [0, A(t)] via E{t) = A{t) in temporal gauge. 
Using the representation of the 7^ in terms of the usual 
Pauli matrices crx,y,z-, we get the Hamiltonian form 

idti) = {-'io'xdx + qA{t)eTx + ma„) • ■0 . (3) 

After a spatial Fourier transform, this gives 

idYk = {[k + qA{t)]cTx + mar„) • 0^ = fffc • 0^ , (4) 

where Hk is a self-adjoint matrix with the eigenvalues 

Hk ■ = ±1/TO^ + [fc + gA(t)]2 = ±ilk{t) , (5) 

and real (time-dependent) eigenvectors u^{t) with the 
properties (m^)^ = 1 and u'^ • = 0. 

Now we expand the spinor solution 0fe(f) of Eq. Q 
into these instantaneous eigenvectors 

0fc(i) = ak{t)e~"‘^At)y^+ ^ (5) 

where we have factored out the rapidly oscillating phase 

t 

Fk(t) = J dt' flkit') ■ (7) 

^0 

Inserting the ansatz ([^ into the Dirac equation Q and 
using • u'^ = ■ uff = 0, we find the evolution 

equations for the Bogoliubov coefficients 

d,.(t) = +Sfc(t)/3fe(t) 6+2*^'= W, 

k{t) = -^k{t)ak{t)e-^^^At) ^ (8) 
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where we have introduced the abbreviation 

. . ^ ^ mqA(t) , , 

-fc(t) = Mfc ■ Mfc = -Mfc • Mfc = ■ (9) 

Note that jafep + |/3fep is conserved and can be set to 
unity |afcp + |/3fcP = 1 by appropriate initial conditions. 
Introducing the ratio Rk(t) = ak(t)//3k(t), we can com¬ 
bine the two equations ([^ into one equality 

Rk{t) = Ek{t) , (10) 

which is the Riccati equation. In contrast to the bosonic 
case (where — |/3/cp = 1), we have a plus sign within 
the brackets, which reflects the fermionic nature of the 
Dirac equation (|afep -I- |/3fcp = 1). 

III. WKB ANALYSIS 


Assuming that A{t) and thus Sfc(t) as well as ipkit) are 
analytic functions in a strip of the complex plane includ¬ 
ing the real axis, the above integral can be calculated 
by deforming the integration contour into the complex t 
plane. In the upper complex half-plane Im(t) > 0, the ex¬ 
ponential decays rapidly—which corresponds to 

the fact that R’^^ is exponentially suppressed. However, 
we cannot deform the integration contour arbitrarily far 
since there will be singularities, e.g., at 


= 0 k + qA(tX) = ±im. 


(14) 


The solution of Eq. (121 can thus be estimated by 


~ exp {2iipk{tl)} ~ exp{-2Im[v?fc(t;()]} . (15) 


Let us consider the most simple example: a constant elec¬ 
tric held E = const. In this case, we have A{t) = Et and 
thus we only get two singularities 


In a stationary situation, the Bogoliubov coefficients 
ttfc and /3fc in Eq. ([^ correspond to solutions with posi¬ 
tive and negative energies, respectively. Thus, j3k can be 
interpreted as the amplitude of an electron in the Dirac 
sea while Uk describes an electron in the positive con¬ 
tinuum. Assuming that the held is switched off initially 
and Anally E{t ^ ±oo) = 0, electron-positron pair cre¬ 
ation can be understood as the transition of an electron 
from the Dirac sea into the positive continuum induced 
by the electric held E{t). Thus, with the initial condi¬ 
tions a™ = 0 and /3“ = 1, the pair-creation probability 
(for the mode k) is given by 


pe+e 
j i. 




ii?rf+ 1 


( 11 ) 


Note that this picture is a bit different from the usual 
description of particle creation in the bosonic case, where 
the probability is given by instead. However, for 

fermions, the particle-hole duality (which is the basis for 
the Dirac sea construction) implies that ak and Pk behave 
symmetrically. Thus, one can use both representations 
(Dirac sea filled with either electrons or positrons) with 
the correct initial conditions. 

Since this probability ® is realistically very small, 
we may approximate |q!™*| « iR'k^^l ^ 1- solu¬ 

tions with small amplitude Rk{t) <C 1, the Riccati equa¬ 
tion (10 \ can be approximated by 


Rfc(f)«5fc(t)e2*^'=W. (12) 


. ±im — k 


(16) 


Insertion into the estimate then yields the usual 
Schwinger exponent in Eq. (IT I. 

In the more general case with more singularities 
they all contribute to the estimate (151 but typically the 
one with the smallest | Im[y:fc(f^)]| dominates (for a more 
detailed discussion, cf. Appendix^. It can happen, how¬ 
ever, that two (or more) have comparable magnitude—in 
which case there can be interference effects, which have 
already been observed in Refs. m UMll El [35], for ex¬ 
ample. 


IV. SAUTER PROFILE 


Let us first apply the method sketched above to the 
case of a Sauter pulse [T] considered previously in Refs. (TJ 
ini. More specifically, we consider an electric held of the 
form 


E{t) =Ei + 


E 2 

cosh^(a;2t) ’ 


(17) 


where E 2 Ei ^ Es and UI 2 m. The associated 
vector potential reads 


Ait) — Eit -\- 


E 2 

— tanh(a;2t), 
0 J 2 


(18) 


Note that this approximation does in general not yield 
the exact prefactor m, but it does give the correct 
exponent—which is the quantity we are interested in. 

Equation (121 can now be integrated with the initial 
condition i?™ = i7^(t —>• — 00 ) = 0 


+ 00 


DOUt 

JrCu 


= Rk{t^oo)^ f df5fc(t)e2*^'=(*). (13) 


and has poles at uj 2 t = i7r/2 -|- iirZ. In terms of the 
dimensionless quantities r = uj 2 t, £ = E 2 /E 1 , p = k/m, 
as well as the combined Keldysh parameter [7J 155] 


muj2 


the equation (141 for the singularities reads 


(19) 


— 00 


T* + e tanh r* = 7(±i — p ). 


( 20 ) 
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Consequently, for small e, we basically get the same sin¬ 
gularity (161 in the upper half-plane as for a constant 


field, which we call the regular singularity 
= 7(* -P)+ 0{e). 


( 21 ) 


However, from the analytic continuation of the vector 
potential ( 181 , we get additional singularities from the 


poles of the complex tanh function, for example at 


^add 2 ^ ' 


( 22 ) 


Now, setting p = 0 for simplicity, we have the follow¬ 
ing picture: For small Keldysh parameters [Eq. @ 1 , 
the regular singularity is closer to the real axis and thus 
its contribution © dominates. In this case, we obtain 
basically the same exponent as for the constant strong 
field El alone (ordinary Sauter-Schwinger effect). If the 
Keldysh parameter exceeds the critical (threshold) value 
of 


Tcrit — 


(23) 


however, the additional singularity = in 12 caused by 
the Sauter pulse E 2 / cosh^ (w 2 i) is closer to the real axis 
and thus its contribution starts to dominate. Thus, the 


exponent in Eq. (151 is reduced and hence the probability 


enhanced (dynamically assisted Sauter-Schwinger effect). 
Inserting the singularity r* into the estimate (151 yields 


Eq 


— +P 
7 


with the auxiliary function 13 [T7] 

(j){z) = z\/1 + -\- arsinh 2 :. 


(24) 


(25) 


For 7 < 7 crit, the regular singularity r*eg « 7(1 - p) 
dominates and we get = f 7 r /2 which reproduces 

Eq. 0. For 7 > 7crit) on the other hand, we have to 
insert = in/2 which yields the results in Ref. [3- 
For Keldysh parameters near the critical value ( [2^ the 
competition between regular and additional singularities 
leads to p-dependent interference effects [HI 113. 


V. GAUSS PULSE 

Now let us repeat the same analysis for another bell¬ 
shaped curve which is visually almost indistinguishable 
from the Sauter profile 1/cosh^(w 2 t): a Gauss pulse 

E{t) = Ei+E2eyip{-{uj2tY} , (26) 

with the same hierarchy as before, i.e., E 2 Ei ^ Es 
and UJ 2 rn (see Ref. [T^ for some numeric results). 
The vector potential can be expressed in terms of the 
error function 

A{t) = Eit + erf(cu 2 t). (27) 

2ui2 


As an important difference to the Sauter profile, the 
Gauss pulse and the error function are analytic in the 
entire complex plane and thus do not give rise to addi¬ 
tional singularities. Hence the only singularities are given 
by Eq. pTl 


T* + £ erf T* = 7 (±i — p) ■ 


(28) 


Again, for small and moderate values of 7 , we basically 
obtain the same regular singularity as in Eqs. (16 \ and 


(211 and thus the same exponent as for the constant 


strong field Ei alone (ordinary Sauter-Schwinger effect). 
Moreover, there appears an infinite set of additional sin¬ 
gularities further away from the real axis whose contri¬ 
butions are therefore suppressed. 

If 7 becomes large enough, however, the regular sin¬ 
gularity starts to deviate from its ordinary position (211 


while its imaginary part approaches that of the additional 
singularities (for more details, cf. Appendix [AT]). Setting 
p = 0 for simplicity, the asymptotic behavior of the error 
function for large imaginary arguments 133 . 


erf ( 17 ) 


exp{7^} 

7V^ 


(29) 


implies that this motion begins when exp{ 7 ^} is large 
enough to compensate the smallness of e, i.e., the critical 
value of the Keldysh parameter roughly scales as 


Tcrit 




(30) 


As a result, although the Gauss pulse looks very similar 
to the Sauter pulse, they behave remarkably differently 
regarding the dynamically assisted Sauter-Schwinger ef¬ 
fect. In one case (Sauter pulse), the critical frequency 
where the enhancement sets in is basically independent 
of £ = E 2 /El whereas in the other case (Gaussian), this 
threshold frequency scales as ~ i/] lne|. 

Apart from the motion of the singularities, the addi¬ 
tional Gaussian pulse does also affect the phase function 
ifk in Eqs. 0 and ( [T^ . Again setting p = 0 for simplic¬ 
ity, we get to first order in £ 




t^Es 

4^1 


1 — ee^ [ It 


-E 


(31) 

where Ii, are the modified Bessel functions of the first 
kind. Note, however, that this expansion is only reliable 
for Keldysh parameters well below the critical value (301. 


VI. OSCILLATION 

As our third example, let us consider a harmonically 
oscillating weak field (see also Refs. | 12 [ [T51 [HI (HI \T7\ 

m) 


E{t) = El + E 2 cos{uj2t ), (32) 
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again with E 2 Ei Es and uj 2 m. Similar to the 
Ganssian case, the vector potential 


E 

A{t) = Eit H- - sin(a;2t) 

UJ2 


(33) 


is analytic in the entire complex plane, 
yields the singnlarities 


Again Eq. (14) 


r* + e sinr* = 7(±i — p). (34) 


The solntions of this transcendental equation can be de¬ 
termined graphically (cf. Appendix A 21. If we again set 
p = 0 for simplicity, the dominant, regular singularity is 
purely imaginary and given by 


Im T* + e sinh(Im r*) = 7 . 


(35) 


As in the Gaussian case, small and moderate values of 7 
basically do not affect the singularity (since e 1) and 
the dynamically assisted Sauter-Schwinger effect requires 
7 to be large enough such that the magnitude of sinh7 
compensates the smallness of e. In this limit, we may 
approximate sinh7 « e'''/2 and thus the critical value of 
the Keldysh parameter roughly scales as 


The worldline instanton x^{u) is a closed loop in Eu¬ 
clidean space-time satisfying the boundary conditions 
a;^(0) = x^(l) and the instanton equations of motion 

mx^ = iqaE^yXu , (41) 

with x^ = dx^/du and x^ = d'^x^/du^. Since E^i, is 
antisymmetric, the absolute value of the “four-velocity” 

Xfj,x^ = Xi + xl + xl + xl = , (42) 

is conserved. 


A. Integral form of the instanton action 

In certain cases, for example a purely time-dependent 
electric field E — E{t)ez, the instanton action A can be 
computed in an especially simple way. Let us consider a 
Euclidean four-potential of the form 

iA3{x4) = — f{x 4 Uj ), (43) 

LJ 


lent ^ I lne| . 


(36) 


Again, we may Taylor expand the phase function ipk in 
Eqs. ([^ and (151 for Keldysh parameters well below this 
critical value 


Im [pfc(ffc)] 


T!'Es 

AEi 


1 — 2e cos(p7) 


-^1(7) 


■ (37) 


VII. WORLDLINE INSTANTON METHOD 


where f{z) is a dimensionless odd analytic function such 
that the electric field E = f'{iujt)ez is real. Inserting the 
vector potential (431, the only nonvanishing components 


of are 1^43 = —1^34 = Ef { 1 x 1 x 4 ) and the instanton 


equations (41) reduce to 


X 3 = -'^^f{uJX4)X4 , 


X4 = -I- 


m 
qEa , 


f{uJX4)x3 , 


(44a) 

(44b) 


It is illuminating to address the above scenarios with 
another technique, the worldline instanton method |34l 
1551145] . In contrast to the WKB approach, the worldline 
instanton technique does not yield the momentum depen¬ 
dence (it only provides the total pair-creation probabil¬ 
ity), but it has the advantage that some calculations are 
easier to do. Again, let us begin with a brief review of 
this method. The total pair-creation probability Pe+e- 
is related to the vacuum persistence amplitude (0out|0in) 


as well as ail = X2 = 0. Thus, in order to have closed 
loops a:^(0) = a;n(l), we must have xi = X 2 = 0 and 
hence Eq. ( [4^ simplifies to 

= xl + x^. (44c) 


Equation (44a I can be integrated immediately to give 


X3 = - - f{u}X4) : 
7 


(45) 


T^e+e- =1-I(0out|0i„)|" . (38) 


If we now express the vacuum persistence amplitude 
(0out|0in) in terms of the path integral and perform a se¬ 
ries of transformations and approximations, we find that 
Pe+e- can be estimated in analogy to Eq. (15) via 

Pe+e- -- e-^ . (39) 


Here A is the worldline instanton action in the presence 
of the electric field E represented by the Euclidean vector 
potential A^{x'^) after analytic continuation to imaginary 
time x^ = {x,y,zAt) = (xi,a;2,a::3,2:4) 


1 



(40) 


where we have set the integration constant to zero in 
order to get closed loops. Using (44cI, we arrive at 


X 4 = -E\la? — x\ = 1-^/^(u;x4). 


(46) 


Furthermore, for a field of the form (43), the instanton 
action (40) can be simplified as well 

1 


A ma 


+ -— f du f{uJX 4 )x 3 

^ J 


,2 \ 4 : 4 c) m 


ma -/ duxi — / duxi . 


(47) 
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0 

FIG. 1: Example of an instanton trajectory x^{u) with 
the appropriate symmetries. In this case, the full instan¬ 
ton action A is four times the action along the highlighted 
(or any) quarter of the trajectory. 


Since is an odd function and thus f (xiuS) and 

f'^{x 4 Uj) are even functions, the instanton trajectory has 
a fourfold symmetry: 0:3 —>■ —X3 and X4 —X4. As a 


result, the contribution to the action (471 is the same for 
each quarter of the trajectory; see Fig.TTr In each of these 


quarters, X 4 is a unique function [Eq. (461] of X 4 and thus 
we can turn the u integral into an X 4 integral 


X4(1/4) 


A = 4- / duxl 

a I 


4m J dx 4 - ^p{uJX 4 ) 

0 3:4(0) 

= 4^ ] 


(48) 


where we have substituted y = qEx 4 /m in the last 
line. The integral coincides (up to substitutions) with 
the function 3(7^) in Ref. (12]. The turning point y* is 
given implicitly by the condition 


7 = /(7X*) ■ 


(49) 


which corresponds to in Eq. (141 for p = 0, where 
we have t* = iyy*. The expression (481 is generally 
much easier to evaluate numerically—and approximate 
analytically—than solving the equations of motion (411. 


VIII. DYNAMICAL ASSISTANCE 

In order to study the dynamically assisted Sauter- 
Schwinger effect described above, we choose the dimen¬ 
sionless function /(y) from Sec. VII A|to be 


fix) = X + e5(x), 


(50) 



-12 -10 -8 -6 -4 

log(f) 

Gauss profile 



Sauter profile 



FIG. 2: Instanton action A ranging from nEs/E (light 
yellow) to 2Es/E (dark blue) for the considered profiles 
as a function of 7 and e. This thin (black) lines are 
contour lines and the additional thick (red) lines depict 
the corresponding estimates (571, (54), and (23), respec¬ 


tively. Note the different scales of the 7 axes. 
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with e ^ 1 and a function g{x) representing the addi¬ 
tional weak time-dependent profile 


E = E[1 + eg'{iujt)] ( 


(51) 


Assuming that g{x) is a well-behaved function (which is 
the case for the Gaussian and the oscillating profile—but 
not for the Sauter pulse due to its singularities), we can 
expand the integral (48 1 to first order in e 




\/l- X + e 


. giix) 



(52) 


B. Oscillating profile 

In order to represent an oscillation, we choose 
g{x) = sinhx- 

In this case, the G{x) in Eq. ([^ becomes 


G( 7 )= 2/1(7). 


(56) 


(57) 


Again, comparison with Eq. (371 shows agreement for 
p = 0 as well as for p = ±2TFpy and p = ± 471/7 etc. 
These additional p values result from the periodicity of 
the oscillation which implies that the times t and t ± T 
with T = 27:1 uj are equivalent. During one period T, 
the acceleration by the strong electric field E causes a 
momentum shift of Afc = qET — qE27:/ui = 27rm/7, i.e., 
Ap = 277 / 7 . 

In view of the asymptotic behavior G{x) 
for large 7 , the threshold for the dynamically assisted 
Sauter-Schwinger effect is reached when 


^3/2g-7 




(58) 


To first order in £, the dependence of x* ^ does not 
contribute because the integrand vanishes at that point. 
The contribution to zeroth order in e reproduces the 
constant field in Eq. (0. The form of G(x) then de¬ 
termines how much, for an additional pulse sg'{iu}t), a 
given strength e and Keldysh parameter 7 , the instanton 
action A is reduced and thus pair production improved. 


A. Gaussian pulse 


For a Gaussian pulse, the Euclidean vector potential is 
represented by the imaginary error function 


gix) = ^ erfix. 


and thus the function G{x) in Eq. (52) reads 

G ( 7 ) = ’ Je ^^/2 ( 




-E 




(53) 


(54) 


Gomparison with Eq. (311 shows agreement for p = 0, 


which demonstrates that the pairs are produced predom¬ 
inantly with p = 0 (as expected). Of course, this is no 
accidental coincidence as the manipulations in Eq. (52) 


for 


are equivalent to those leading to Eqs. (311 and 
p = 0; see also Ref. pS] • 

For large 7 , the function G(x) behaves as expjy^j/y^ 
and thus the threshold for the dynamically assisted 
Sauter-Schwinger effect is reached when 


3 -A 

7 g r 


(55) 


which is consistent with Eq. (30l. Again, we would like 


to stress that the first-order expansion (541 is only valid 


for Keldysh parameters 7 well below this critical value. 


which is again consistent with Eq. (36 1 . 


IX. STANDING WAVE 


One might perhaps object that the fields E{t) con¬ 
sidered here are not solutions of the Maxwell equations 
in vacuum, i.e., without sources. This objection could be 
motivated by the idea that the fermionic modes somehow 
“feel” the sources generating the electromagnetic back¬ 
ground field and thus behave totally different in field 
configurations with and without sources. However, re¬ 
tracing the steps of the formalism presented above, we 
see that we only used the Dirac equation in a given ex¬ 
ternal background field A^. Whether this field is a 
solution of the Maxwell or Proca or Yang-Mills equations 
with or without sources is irrelevant for our derivation. 
Note that we also did not exploit local gauge invariance. 

In order to illuminate this point, let us consider the 
following electric field profile: 

E{t, x) = [El ± E 2 cos{ujt) cos(A:a:)] . (59) 


It can be represented by the Euclidean vector potential 


iAzixi^Xi) 


E 1 X 4 + 


E 2 

bJ 


sinh(a;a; 4 ) cos(fca;i), 


(60) 


which generates the above electric field in z direction as 
well as a magnetic field in the y direction. Thus, the only 
nonvanishing components of E^^, are 


7^43 = —iF 34 = El + E 2 cosh(a;a; 4 ) cos(fca;i), ( 6 Ia) 

k 

iFi 3 = —iF^i = —E 2 — sinh(a;a: 4 ) sin(fca:i), (61b) 
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which yields the following instanton equations: 


mil = — qaE2 — sinh(wa:4) sin(fcxi)i:3 , ( 62 a) 

LJ 

mx 3 = — qa [Ei + E 2 cosh(a;x 4 ) 003 (^ 0 : 3 )] ±4 
k 

+ qaE2 — sinh(wa:4) sin(fca;)ii, ( 62 b) 

OJ 

mxi =qa [Ei + E2 cosh(a;a;4) cos(fca;3)] ±3 . ( 62 c) 


The equation (62a I for the transverse component xi is 
solved by xi(u) = 0 or a:i('u) = ± 7 r/A: etc. Choosing 
a solution in the maximum of the electric field such as 
xi (u) = 0 , the other equations reduce to 

mis = —90 + E 2 cosh(a;a; 4 )] X4 , (63a) 

mx 4 = qa [Ei + E 2 cosh(a;a: 4 )] is . (63b) 

As a result, we find that the instanton trajectories and 
thus also the instanton action A are not affected by the 
additional factor cos{kx) in Eq. (59) at all. This invari¬ 
ance can be generalized to vector potentials ^ 3 ( 0 : 1 , 0 : 4 ) 
for which dA^/dxi vanishes at xi = 0 for all X4, which is 
the case for all even functions A^{—xitX 4 ) = ^ 3 ( 0 : 1 , 0 : 4 ), 
for example. 

Previous experience suggests that temporal variations 
of the electric field E{t) tend to increase the pair-crea¬ 
tion probability whereas longitudinal spatial variations 
E{x)ex tend to decrease it (an interplay between these 
two tendencies is discussed in Ref. m-, see also Ref. |33] ). 
Our example shows that a transversal spatial dependence 
does not necessarily have this effect because the worldline 
instanton action A is actually independent of k. Note, 
however, that the situation would change drastically for 
a propagating wave such as cos(wt — kx) instead of a 
standing wave cos(wt) cos(fca;); see, e.g.. Ref. jTU]. Fur¬ 
thermore, the prefactor in front of the exponential in 


Eq. (39), which can be derived via studying small per¬ 


turbations around the instanton trajectory, could well 
depend on k since it affects the effective pair-creation 
volume. 


For uj = k, the field (59) represents a vacuum solution 


of the Maxwell equations in the form of a standing wave, 
which can be generated by the superposition of two plane 
waves cos{ujt ± kx) with the same polarization but prop¬ 
agating into opposite directions. For uj ^ k, however, 
this field profile only solves the Maxwell equations with 
nonzero sources. Nevertheless, the instanton action A 
is the same for uj = k and lo ^ k, which demonstrates 
that it is irrelevant whether the background field A^ is 
a vacuum solution of the Maxwell equations or not. The 
difference between a standing cos(ujt) cos{kx) and a prop¬ 
agating wave cos(ujt — kx) is far more important than the 
question of whether the field is a vacuum solution of the 
Maxwell equations (w = k) or not {uj ^ k). 

Of course, we do not deny that it is certainly desirable 
to calculate the pair-creation probability for field config¬ 
urations which are as realistic as possible. On the other 
hand, our results show that one should not discard results 


just because the considered field profile is not a vacuum 
solution of the Maxwell equations m- 


X. CONCLUSIONS 

Via the WKB and the worldline instanton method, we 
studied and compared the dynamically assisted Sauter- 
Schwinger effect for three profiles (Sauter, Gauss, and 
oscillating) and found qualitative differences. For the 
Sauter pulse considered previously, the critical (thresh¬ 
old) Keldysh parameter (23) is determined by the compe¬ 
tition between the regular singularity (211 and the addi¬ 


tional singularity (221 and is basically independent of the 


magnitude E 2 of the weak field. By contrast, the other 
two cases do not feature such a competition between dif¬ 
ferent singularities and the critical (threshold) Keldysh 
parameters (301 and (36 1 do depend on the magnitude 


of the weak field, albeit only logarithmically. These pro¬ 
found differences can be traced back to the distinct fea¬ 
tures of the three profiles in the complex t plane—which 
also determines the momentum dependence including in¬ 
terference effects etc. 

Furthermore, one would expect that these differences 
do also affect the behavior of the prefactor (which was not 
considered here) in front of the exponential in Eqs. (15) 


and (39), e.g., its dependence on e. Note, however, that 


this prefactor cannot resolve or counterbalance the ob¬ 
served difference between the pulse profiles (e.g., Sauter 
and Gauss). For example, let us consider a Keldysh pa¬ 
rameter a bit above the threshold value ( [ 2 ^ of 7 crit = 
tt/2. For a Sauter pulse, the dynamically assisted Sauter- 
Schwinger effect is already operational in this regime— 
which has also been confirmed by numerical computa¬ 
tions of the total pair-creation amplitude (including the 
prefactor); see, e.g.. Ref. [TB]. For a Gaussian profile with 
the same parameters, on the other hand, the weak-field 
contribution (x E 2 is negligible in the vicinity of the in¬ 
stanton trajectory for small values of e. Since the prefac¬ 
tor can be determined by considering small perturbations 
around the instanton trajectory, it is thus also almost 
unaffected by this weak-field contribution cx E 2 . This 
behavior was confirmed by preliminary numerical inves¬ 
tigations [32] which indicate that the Gaussian prefactor 
is approximately constant (i.e., independent of e and 7) 
for small 7 < 1 but then starts to grow with both, in¬ 
creasing e and 7, for larger values of 7. By contrast, the 
prefactor of the Sauter pulse is nearly independent of e 
for small as well as large 7 and only varies with e near 
threshold 7 « 7crit = 7 r/ 2 . (Note that one cannot simply 
take the limit e —^ 0 here since the instanton approx¬ 
imation would break down eventually.) In both cases, 
however, the moderate variation of the prefactor is sub¬ 
dominant in comparison to the strong enhancement of 
the pair-creation probability induced by the changes in 
the exponent considered here. 

Our results can be generalized in a straightforward 
manner to other profiles, for example a Lorentzian pro- 













file, which is expected to behave similarly to the Sauter 
pulse. The fact that two visually almost indistinguishable 
bell-shaped pulses (Sauter and Gaussian) display such 
drastic differences demonstrates that tunneling in time- 
dependent backgrounds can show surprising and partly 
counterintuitive results—which motivates further studies 
in order to understand these nonperturbative phenomena 
better. 
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Appendix A: Positions of the singularities 


The outgoing ratio R™*, which determines the pair- 
creation probability, is approximated by the contour in¬ 
tegral in Eq. according to the WKB method. In 
order to find a suitable integration contour in the upper 
complex half-plane Im(t) > 0 , we need to know the posi¬ 
tions of the integrand’s singularities given by Eq. (141. 
These singularities are poles of the function 5^ (t) as well 
as square-root-type branch points (with associated cuts) 
of (pk{t). The weak electric field may introduce further 
complex singularities into the vector potential A{t) itself 
and thus into the integrand; however, this is neither the 
case for the Gauss pulse nor for the oscillation. 


these additional singularities. In this limit, the r*- 


independent terms on the right-hand side of Eq. (28 1 can 


be neglected, and we find the series of singularities 


^add,r; 


Ine 




cos( 20 „) 


(A2) 


with n G Z. The complex phase is the single real 
solution of the transcendental equation 


I lne| tan(26>„) — 20 „ = 27rn 


(A3) 


in the range 7 r /4 < < 37 r/ 4 , which can be found nu¬ 

merically. Note that, in this approximation, the depen¬ 
dence on momentum p and the difference between the 
plus and the minus cases in Eq. (28) disappear. We 


thus only get one result where there are actually two 
different singularities (± cases). For fixed e, all points 
r*jjd n ^ hyperbola in the complex plane since 

Im^^T-add,™) - Re^('ra*dd,n) = Une| for all n. If we con¬ 
sider a fixed n instead, we can approximate the solution 


of Eq. (A3) for sufficiently large | lne| by linearizing the 


tangent around 29 ^ = tt. To leading order, this gives the 
asymptotic form 


‘ add,n 


7 r(n -I- 1/2) 
\/|lne| 


f\/| ln£ 


(A4) 


which parametrizes another hyperbola as e is varied, cf. 

Fig.M 


For p = 0, the singularity equation (28) is symmetric 
with respect to the imaginary axis (Rer* —)■ — Rer*) 
and the regular singularity remains on this axis for all 
values of e and 7 . The movement of the singularities is 
depicted in Fig. 


1. Gauss pulse 


2. Oscillation 


The dimensionless singularity positions t* = 012 are 
given by Eq. ( p^ in this case. The regular solution 
[cf. Eq. (HH, which is the only solution to stay finite in 
the limit £ —)■ 0, can be Taylor expanded around £ = 0 


R*eg = 7(*-p) - ^erf[7(i-p)]£-k 


X erf[ 7 (i — p)]£^-I- C>(£^). (Al) 


For values of 7 below the threshold 7 crit ~ \J\ ln£|, the 
terms of order 0{e) and higher are much smaller than 
unity. This is no longer true, however, as 7 approaches 
7 crit- Hence, the truncated Taylor series (Al) is appropri¬ 


ate to describe the movement of the regular singularity as 
the parameters are varied within the subcritical regime. 
The other solutions of the singularity equation (28) 


diverge as e approaches zero, so we may approximate 
the error function by its asymptotic form erf r ^ 
— exp(T^)/(y^r) [37] for small £ in order to describe 


a. Graphical solution of the singularity equation 


The transcendental equation (34) for the dimensionless 


singularities r* can be solved graphically. We start to 
develop the graphical solution scheme by splitting the 
equation and the complex variable t* = u* + iv* into 
real and imaginary parts 


u* -I- £ sin(u*) cosh(u*) = —pj , 
V* + e cos{u*) sinh(u*) = ±7. 


(A5) 

(A 6 ) 


For singularities in the upper complex half-plane (u* > 0) 
with sin u* ^ 0, Eq. (A5) can be solved uniquely for the 
imaginary part v*, which yields v* = Using (u*) with 

Using(u) = arcosh(^— ^ +P7 \ ^ 

\ £smu/ 


Inserting this expression into Eq. (A 6 ) leads to 
F{u*) = ±1 


(AS) 
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Im T* 



(a) Singularities moving as 7 is varied from 1 to 7 (dark red to light 
yellow dots) for constant e = 10“^. The associated threshold for 
the dynamic assistance is yj\ In e\ 4. In the subcritical 7 range, 
the regular singularity on the imaginary axis moves essentially as 
in the case of the ordinary Sauter-Schwinger effect (r*gg 5 ^ ^ 7 ), 
cf. Eq. ( |16| l, and the additional singularities are much deeper in 
the complex plane. As 7 gets close to the critical threshold, the 
regular singularity slows down, cf. Eq. and “pushes” the 

additional singularities away from the real axis. 


Im T* 



(b) Singularities for 7 = 4 and e = 10® with s going from —11 
to —3 (dark red to light yellow dots). The black crosses are the 
corresponding asymptotic solutions ^ for n = — 1,..., —4, see 
Eq. ( |A2[ |, which move along the dashed hyperbolas as e is varied. 
The transition to the critical regime is at £ 5 ^ 10“^. At this value, 
the regular singularity on the imaginary axis starts to move to¬ 
wards the real axis. The additional singularities approach the real 
axis continuously as e grows. In the subcritical range, the asymp¬ 
totic values ^ are appropriate to describe the movement of 
the additional singularities. 


FIG. 3: Numerical solutions r* of the Gauss-case singu¬ 
larity equation (281 for the momentum p = 0. 


with the function 


F{u) = - 


arcosh — 


M + P7 

£sinM 


cos{u)] 


u + p'j 
sin It 


— 


(A9) 


F{u) 




(b) Imaginary parts of the singularities over their real parts. 


FIG. 4: Graphical solution of the singularity equa¬ 
tion (341 for p = 0 and £ = 10“^, so 7crit is roughly 
|lne| « 11.5, cf. Eq. (36l. The function F{u) is even 
for p = 0 , so we concentrate on positive u for simplicity. 
The real parts u* of the singularities r* are determined 
in the upper plot 0 . For p = 0, F(tt > 0) is only real 
for u G (tt, 27r), (37r,47r), etc. and increases strictly over 
each one of these intervals, respectively, so there are two 
singularities per interval. The corresponding imaginary 
parts V* are found by evaluating Using (u) shown in the 
lower plot (0 at M*. 


The real equation (A 81 can be solved graphically (see 
Fig.0 in order to find all singularities r* with real parts 
satisfying sinu* ^ 0. The corresponding imaginary parts 
of the singularities are given by Using (u*). 

Let us now treat the special case where sinu* = 0, 
i.e., those singularities with real parts u* G ttZ, which 
are not covered by the graphical solution scheme. Then, 
Eq. (A51 reduces to u* = —P 7 , so this case can only 
occur for particular momenta p G 'L'k The cosine in 
the lower Eq. (A 6 ) is cosu* = ±1 for these singulari¬ 
ties. For cosu* = 1 {u* = 0, ±27r, etc.), Eq. (A 6 l has 
one unique solution u* > 0 fulfilling u* -I- ssinhu* = 7 . 
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But for cosu* = —1, in contrast, the left-hand side of 
Eq. (A 6 l does not increase strictly with respect to v*, so 
this equation can have multiple solutions. It turns out 
that there are three different positive solutions v* for (ap¬ 
proximately) subcritical 7 < 7 crit ~ |lne|, see Eq. (36 1 , 
and one positive solution for supercritical 7 . In the case 
of multiple singularities r* with the same real value u*, 
the singularity which lies closest to the real axis domi¬ 
nates. 


b. Movement and asymptotics of the singularities 


In the subcritical regime, the regular singularity is well 
described by the Taylor series 

£2 

A*eg = 7 (*-P) -sin[ 7 (i-_p)]e-Fsin[ 27 (i -p)]—+ 0 {e^). 

(AlO) 

As in the Gauss pulse case, the singularity equation ( [M| 
can be solved asymptotically for p = 0 and e —)• 0 (thus 
Imri^jd 00 ) in order to describe the additional singu¬ 
larities. Here, we get 


add.i 


7r(2n — 1) -I- i| lne| 


(All) 


with n G Z, i.e., each additional singularity converges 
against an edge of its “tt interval” for small e\ see Fig. 
Note that each n represents two different solutions of 
the singularity equation (± cases) again since these two 
solutions merge asymptotically. 

If we want to sum up the contributions from all singu¬ 
larities to the contour integral (131, we have to approx¬ 


imate the singularities for fixed parameters p, e, and 7 
instead of fixed n. Concentrating on the mode p = 0 
for simplicity and inspired by the graphical method in 
Fig. a we find that the real parts u* of the additional 
singularities do asymptotically (|u*| 3> 1 ) coincide with 
the centers of the “tt intervals”, i.e., T*dd ~ '’’i « for large 
n S N with 


Im T* 



(a) Singularities for constant e = 10“® (7crit ~ 11.5) and 7 
growing from 7 to 20 (dark red to light yellow dots). The lin¬ 
ear movement of the regular singularity Ri i-y (ordinary 

Sauter-Schwinger effect) slows dow n exponentially near the criti¬ 
cal threshold according to Eq. | |35[ l. The additional singularities 
move along the complex curves u-|-iuaing(u), see Fig. |4(b)| which 
are independent of 7 for p = 0 . 

Im T* 


25 

20 

15 

10 

5 

0 



n 2 n 




Re T* 


(b) Moving singularities for 7 = 13.8 and e = 10® with s = 
— 10,..., —2 (dark red to light yellow dots). The critical thresh¬ 
old is at £ Ri 10“®. The regular singularity starts to approach the 
real axis as the nearby additional singularities come close. For de¬ 
creasing £, each additional singularity approaches the ima ginary 
axis, converging towards the edge of its “tt interval”, cf. Eq. ( |aIT| i. 


T-i.n = ±7r ( 2n- - 


= ±7r I 2n - - 


■ losing (±27rn=F tt 


/27rn — n/2\ . . , 

i arcosh ( -j . (AI2) 


In this context, the ± sign represents singularities with 
positive/negative real parts. Apart from that distinction, 
each n corresponds to two actual singularities again: the 
± cases of Eq. (A 8 I, which merge in the considered limit. 
Since arcosh a; » ln(2a;) for a: ^ 1, the imaginary part of 
the additional singularities grows logarithmically with n 
while their real part increases linearly. 


Appendix B: Calculation of the contour integral (13) 


In order to estimate the integral (131, the integra¬ 


tion contour can be shifted away from the real axis, up¬ 
wards into the complex plane, until the first singularity 


FIG. 5: Numerical solutions t* of the singularity equa¬ 
tion (341 for the oscillating weak field and p = 0. 


is met at some Since the integrand is analytic for 
0 < Imt < Imt^, this operation is permitted and the 
behavior of the complex phase ipk(t) gives an estimate 
of the integrand Sfe(f) exp[ 2 zy)fc(t)] with now complex t. 
However, for the estimation of the integral, additional 
control of the integration along the line Imt = ImtJ is 
required, which is rather difficult to obtain. Therefore, 
the integration contour is shifted further upwards, be¬ 
hind the singularities (see Fig. |^, and Gauchy’s theorem 
is applied. 

The phase ipk{t) and the integrand can be analyti¬ 
cally continued, however, locally on different Riemann 
sheets (with cuts in between) as the singularities are of 
a square-root type. The contribution from such a shifted 
integration contour is subdominant due to the exponen- 
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Im t 



FIG. 6: Example for the deformed complex integra¬ 
tion contour (dotted parts: no contribution; solid parts: 
branch cut contributions and residua) in the case of three 
singularities ]^_3 (crosses) with branch cuts running ver¬ 
tically upwards (dashed lines). 


tial damping coming from the phase (pk{t) [Eq. 01 , but 
there remain contributions from integration along both 
sides of the cuts and around the singularities. In order to 
understand how they can be controlled, it is instructive 
to consider first the ordinary Sauter-Schwinger effect as 
an example. 


c. Contour integral for a constant electric field E 


The advantage of this example is that here contribu¬ 
tions from the singularity and from the cut are explic¬ 
itly determinable. Since the singularity Vf in Clkit )—the 
plus case of Eq. (161 is relevant in this context—is of 
square-root type, the phase which is obtained as 

an integral of Clk{t), behaves as (t —The residuum 
(representing a circle integral around the singularity) can 
be calculated exactly 


Sp = —— exp 


Es 

1^ 


(Bl) 


with the function (l){z) defined in Eq. (25l. The complex 


cut emerging from that singularity, chosen vertically up¬ 
wards from produces a jump in the phase function 
ipk{t), however, only in its real part. In consequence, the 
absolute value of the integrand Sfc(t) exp[2ji^fc(t)] is con¬ 
tinuous across the cut and decays merely polynomially 
for Imt > Imf^ due to 5fc(t); it does not decay expo¬ 
nentially along the cut as might be expected from the 
behavior for Imt < Imt^. However, the difference of 
values on both sides of the cut oscillates rapidly. The 
contribution from the cut can be written as 


Sc = exp 


(f+i«p))§ 


Ic{u) du (B2) 


J Jo 


where Ic{u) « s\n[2u^Es/{^E)]/u is a fast-oscillating 
function of u. The integral converges solely due to its 


fast oscillation and gives approximately the number tt/6. 
In effect, both contributions to 


i?™* = Sp + Sc ^ - ^ exp 


Es 

E 


, (B3) 


from the singularity and along the cut are of the same 
order. We take it as suggestion that neglecting the 
contribution from the cuts in more complicated (time- 
dependent) electric fields should not change the final re¬ 
sult by more than a constant numerical factor. 


d. Contour integral for the oscillating weak field 


In the case of the oscillating-weak-field profile given 
by Eqs. (32 h( 33), there are infinitely many singularities 

Considering p = 0 again. 


as shown in Appendix |_ 

we get two additional singularities per “tt interval” (see 
Fig. the singularities with Re < 0 are just mirror 
images of those with positive real parts) plus the regular 
singularity on the imaginary axis. In order to calculate 
the contour integral we need to integrate around all 
these singularities (generating residua) and along both 
sides of the corresponding branch cuts. As the example 
above shows, it can be expected that the contributions 
from the cuts to i?™* are of the same order as the residua 
S'p(f^) at the singularities = t* jooi. Hence, i?™* is 
approximately given by the residual sum 


DOUt 


Spir* ) + V 5'p(r*dd) 


(B4) 


in dimensionless time units t = uj 2 t. At each singularity, 
5fc(t) has a pole while pk{t) has merely a square-root- 
type branch point, so pk{t) may be treated as constant 
when integrating 5fc(f) exp[2i(^fc(t)] along a small circle 
around the pole. This way, we get 

Sp{T*) = ( 35 ) 


where the sign depends on whether r* is a plus or mi¬ 
nus solution of the singularity equation (34) and on the 
concrete potential A{t). Although we cannot sum the 
series (B4) exactly, we can prove its convergence and es¬ 
timate its value. Using the asymptotic form of the singu¬ 
larity positions T*dd for I Rer*ddl ^ 1 and p = k/m = 0 
[Eq. (A12)], we can approximate the values of po at most 
of the singularities and estimate the sum over all addi¬ 
tional singularities appearing in Eq. (B41 


T* 

add 


<Ei^p«dd)i 


^ ^ g-2Im[¥3o(TUd)] _ 

T* 

add 

(B6) 
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Taking into account the contributions from both series of 
singularities, i.e., the ± cases in Eq. (341, we get 


g-2Im[¥3o«dd)] 


OO 

n—1 

+ exp{-2Im[(/3o('r*,„)]}^ (B7) 


^exp{ - 2 Im[(po «,„)]} 


Inserting it into the above sum and performing some fur¬ 
ther minor simplifications, we find the estimate 


g-2Im[¥3o«dd)] < 

"^add 


4(^) 


TvqEi /a;| 


f SnEi \ 

V J 


ATzqEi j ijj^ 


- 1 


(B9) 


in terms of the usual physical quantities. The right-hand 
side is finite, which proves the convergence of the residual 
sum (at least for p = 0). As a check, two limiting cases 
can be considered: 


with „ 


from Eq. (A12). 


For singularities with large 


real parts |Rer*| 3> 1, we neglect the rv? term under 
the square root in Ofc(t) [see Eq. (§1 , and after some 
other minor analytic simplifications we get 


Im[po(r*)] 


Es 


ReT*|(Im T* 


I). 


(B8) 


a) El —)■ 0: As the assisting oscillation vanishes, we 
expect the ordinary Sauter-Schwinger effect. Con¬ 
sistently, the right-hand side of Eq. (B91 becomes 
zero in this limit—there is no contribution from the 
additional singularities. 

b) W 2 —t 0: If the oscillation becomes a static field 
which has a negligible influence due to E- 2 , ^ Ex, 
the right-hand side of Eq. (B91 vanishes again, leav¬ 
ing the ordinary Sauter-Schwinger effect, as ex¬ 
pected. 
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